Density functional theory for nearest-neighbor exclusion lattice gasses in two and 

three dimensions 



o 
o 

(N 

wo: 
< 



CO 

a 

C 

o 
o 



> 
in 

On 
oo 

o 
m 
O 

c3 



o 
o 



X 
S3 



Luis Lafuent^] and Jose A. Cuestjj] 
Grupo Interdisciplinar de Sistemas Complejos (GISC), 
Departamento de Matemdticas, Universidad Carlos III de Madrid, 
Avda. de la Universidad 30, 28911-Leganes, Madrid, Spain 
(Dated: February 2, 2008) 

To speak about fundamental measure theory obliges to mention dimensional crossover. This fea- 
ture, inherent to the systems themselves, was incorporated in the theory almost from the beginning. 
Although at first it was thought to be a consistency check for the theory, it rapidly became its fun- 
damental pillar, thus becoming the only density functional theory which possesses such a property. 
It is straightforward that dimensional crossover connects, for instance, the parallel hard cube system 
(three-dimensional) with that of squares (two-dimensional) and rods ( one- dimensional) . We show 
here that there are many more connections which can be established in this way. Through them 
we deduce from the functional for parallel hard (hyper)cubes in the simple (hyper)cubic lattice the 
corresponding functionals for the nearest-neighbor exclusion lattice gases in the square, triangular, 
simple cubic, face-centered cubic, and body-centered cubic lattices. As an application, the bulk 
phase diagram for all these systems is obtained. 

PACS numbers: 05.50.+q, 64.60.Cn, 61.20.Gy, 05.20.Jj, 64.10. +h 



I. INTRODUCTION 



Rosenfeld's fundamental measure (FM) theory [lj is 
peculiar among the weighted density approximations 
(WD As), and it is so for many reasons. To begin with, 
classical WDAs are constructed upon the knowledge of 
the thermodynamics and structure of the uniform fluid, 
while FM theory is constructed on a geometrical basis. 
Originally it needed scaled-particle theory to produce a 
functional, but in its latest formulations this is not a 
requirement anymore 0, Q , and geometry stands as its 
unique ingredient. Another important difference is that, 
while the extension to mixtures of classical WDAs is far 
from being straightforward, the natural formulation of 
FM theory is for a mixture (although it has recently been 
shown that a FM theory for hard spheres does not ac- 
commodate a description for mixtures as well as it was 
previously thought One further remarkable feature 
is that FM theory performs best where classical theo- 
ries are poorest: in the high density region. It has been 
shown, for instance, that the description FM theory pro- 
vides for a solid is extraordinarily accurate in all its de- 
tails 0, Q ■ This is probably the reason why the belief 
has spread that FM theory is the best density functional 
theory for the system of hard spheres. But there is no 
free lunch. Such a peculiar structure makes the theory 
extremely rigid, so much that it is very difficult (some- 
times impossible) to improve a particular detail without 
spoiling another. This shows up very clearly if one tries 
to improve the equation of state for the liquid phase. FM 
theory yields the scaled particle equation of state. The 
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large difference in accuracy between the liquid and the 
solid gives rise to a not very good prediction of freezing 
[3(. If one tries to replace the equation of state by, e.g. 
Carnahan-Starling, the internal structure of the theory 
squeaks |6| and looses some of its nice features (although 
for some purposes the defects may be mostly irrelevant 
0)- 

But by far the most genuine property of FM theory is 
dimensional crossover, something that this theory shares 
with the exact functionals and with no other known den- 
sity functional theory. In its first formulation this prop- 
erty was introduced as the ability of a modification of 
Rosenfeld's original functional for hard spheres to recover 
the exact functional for 'zero-dimensional' cavities (cav- 
ities which cannot hold more than one sphere) @, but 
it was immediately extended to describe the property of 
the exact d-dimensional functional to reproduce the ex- 
act (d — l)-dimensional one when evaluated at a density 
profile which is delta- like on a hyperplane ;9|. Needless 
to say that functionals having dimensional crossover are 
particularly suitable for studying fluids under strong con- 
finement. 

The first modified FM functional for hard spheres pro- 
duced accurate functionals for d = 2 and d = 1, apart 
from yielding the exact one for d — (cavities) |9| . When 
applied to the system of parallel hard cubes, whose FM 
functional can be obtained for arbitrary dimension (be- 
ing exact in d = 1) 10j, it was shown that dimensional 
crossover consistently transforms the <i-dimensional func- 
tional into the (d — l)-dimensional one, down to d = 0. 
The acknowledgment that this property lies at the heart 
of the formal theory suggested the last step in this di- 
rection: transforming thisproperty into the constructive 
principle of FM theory [3, 0- Under this new formu- 
lation FM theory has been generalized to systems with 
soft interaction potentials [llj, anisotropic hard-particle 
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models Il2l , nonadditive mixtures , latti ce g ases 
[TH Hil Il7| . and even fluids in porous media 0, 119 ]. 

So we see that dimensional crossover was first looked 
at as a very stringent constraint on density functionals 
and later as a way of rising from d = and d = 1 to 
d > 1 in the construction of FM functionals. But dimen- 
sional crossover has another use which has hardly been 
exploited: one can get new systems out of known ones. 
The first (and to our knowledge the only) example of such 
a use was already provided by Rosenfeld et al. jsj. By 
forcing hard spheres to have their centers of mass on one 
out of two parallel planes separated a distance shorter 
than a sphere diameter, one obtains a binary mixture 
of nonadditive hard spheres with negative nonadditivity 
(2<7i2 < an + &22, with cry the center-to-center contact 
distance between spheres of type i and j). The amount 
of nonadditivity depends on the distance between planes. 

Physically, dimensional crossover amounts to applying 
an infinite-strength external potential all over a (d + 1)- 
dimensional lattice, except in a certain d-dimensional set 
of sites which defines the effective system. But this is not 
the only way to construct functionals for new systems out 
of known ones. One can also apply an appropriate exter- 
nal potential at selected sites and modify the interaction 
accordingly, thus obtaining a new system without reduc- 
ing dimension. This trick has already been applied to 
obtain the exact functional for a nonadditive mixture of 
hard rods in a one-dimensional lattice from that of the 
additive mixture 

One can consider all systems which are related through 
the kind of transformations we have just described above. 
Then, because of the internal consistency FM theory has, 
the functional for a given model "contains" the functional 
for any other model to which it is related. Our goal in 
this paper is to show how this works for a set of well- 
known lattice gases. The methods we will use are very 
general, so their application to other families of models 
should not be difficult. This is interesting if we take into 
account the importance that lattice gases are getting in 
the study of certain inhomogeneous problems, such as 
the behavior of fluids in porous media 123, 12JJ . 

The paper is organized as follows. Section ITT1 describes 
the general procedure to obtain the excess free-energy 
functional for nearest-neighbor exclusion lattice gases in 
different lattices, either starting from a higher dimen- 
sional functional for cubes and using dimensional reduc- 
tion to a plane or a hyperplane (this procedure is sub- 
sequently applied to the square, triangular, simple cu- 
bic and face-centered cubic lattices), or starting from 
the functional of cubes and applying an infinite-strength 
external potential in the appropriate set of lattice sites 
without reducing the effective dimensionality of the sys- 
tem (this procedure is the one applied for the body- 
centered cubic lattice). The final (simple) expressions 
for the functionals are explicitly obtained in closed form. 
In Section lTTTl these functionals are applied to obtain the 
bulk phase diagram for all the systems considered. There, 
FM theory results are compared with those from other 



classical theories, showing that the former is at least at 
the same level of accuracy than the latter. Finally, we 
conclude in Section llVl 



II. THEORY 

In this section, we will derive the FM functional for 
lattice gases with nearest-neighbor exclusion in five dif- 
ferent lattices: square (hard-square lattice gas), triangu- 
lar (hard- hexagon lattice gas), simple cubic (sc), faced- 
centered cubic (fee) and body-centered cubic (bcc). All 
these systems have been already considered in the litera- 
ture as simple models for the hard-sphere system |22J . As 
explained in the Introduction, the derivation of the four 
first will make use of the dimensional crossover property 
of FM functionals, so all the four will be obtained from 
the known functional for (d+l)-dimensional parallel hard 
cubes in a simple (hyper)cubic lattice [l(j,Cj| being the 
dimensionality of the final systems). For the last one, we 
will start from the functional for the three-dimensional 
parallel hard cubes in the simple cubic lattice, and we 
will apply an infinite-strength external potential to the 
appropriate set of lattice sites of the original lattice, so 
that the effective lattice becomes a bcc one. The fact that 
all the models considered exclude only nearest neighbors 
forces the edge length of the cubes to be a — 2 lattice 
spacings. 

Of the different ways in which the FM functional for 
this particular model of (hyper)cubes can be written [Trl 
Il7j . the simplest expression for the excess free energy is 
probably 

mM = Y. E (-1)^*0 (« (k) («)), a) 

sez d ke{o,i} d 
where n( k )(s) are the weighted densities 

n^(s)= £ p(s + r) (2) 

reB(k) 

labeled by the vector index k — (ki, . . . , kd), k = 
J2i = i ki, *8(k) denotes the set 

B(k) = {re {0, l} d : < n < k u i = 1, . . . , d} (3) 

and $o (rj) = rj + (1 — ij) ln(I — ?/) is the excess free energy 
for a zero-dimensional cavity with mean occupancy < 
r/ < 1 (/3 is the reciprocal temperature in Boltzmann's 
constant units). 

A. Square lattice 

The kind of dimensional reduction we have to perform 
in order to obtain the hard-square lattice gas out of cubes 
in a cubic lattice is illustrated in Fig. ^ It amounts to 
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TABLE I: Coefficients of the linear combinations 
Xw°z(k)p(t + e;) appearing in ©, which define the 
weighted densities for the hard-square lattice gas model. 



FIG. 1: Hard cubes with edge length a = 2 in a simple cubic 
lattice. Their centers of mass are constraint to lie on the plane 
si + S2 + 2s3 = 0. Filled circles mark the sites excluded by the 
cubes which lie on this plane. Notice that on the plane, cubes 
behave as the hard-square lattice gas with nearest-neighbor 
exclusion. The diagram at the lower right corner helps to 
visualize the relative position of the cubes in the figure. 



forcing the centers of mass of the cubes to lie in the plane 



as 



E *° ( E S ( Sl + M ) E p(s 2 + r 2 + r 3 ,s 3 + r 3 ) 
-0 (r|u,k) 



st 



where (r\u, k) denotes those vectors r e £>(k) which ver- 
ify n + r 2 + 2r 3 = u. Now, one and only one of the delta 
functions within the sum in r is 1, the others are 0, so 
the above expression admits the following rewriting 



"P sq = {(si,s 2 ,s 3 ) € Z 3 : Sl 



si 



2s 3 = 0}. (4) 



(Notice that we might have chosen several equivalent 
planes, given the symmetry of the system.) Figure ^ 
shows that the effective underlying lattice is a square 
lattice defined, for instance, by the pair of orthogonal 
vectors {ei = b 2 — bi, e 2 = b 3 — bi — b 2 } of the plane 
■p sq , {bi,b 2 ,b3} being the canonical vector basis in Z 3 . 
Furthermore, the effective interaction potential between 
the cubes within this lattice is nearest-neighbor exclu- 
sion. 

In terms of the one-particle distribution function, the 
dimensional reduction can be imposed by setting 



p(s) = p(s 2 + s 3 ,s 3 )8(si + s 2 + 2s 3 ), 



(5) 



where 8(0) = 1 and S(x ^ 0) = (a Kronecker symbol). 
The dependence of the two-dimensional density profile 
reflects the choice of basis vectors. 

To obtain the FM functional for the hard-square lattice 
gas it only remains to introduce the density profile JSJ in 
the functional through the weighted densities J2J an d 
to eliminate the unnecessary degrees of freedom. In what 
follows, we will carry out this task in detail. 

The (signless) contribution of one weighted density, say 
m k )(s), to the excess free-energy functional will be 



E $0 E p(s 2 + r 2 +r 3 ,s 3 +r 3 )8(s 1 +u) 

sSZ 3 \reB(k) 

where u = n + r 2 + 2r 3 and we have made use of the 
translational invariance in s. A better way of expressing 
this is to split the sum in r according to the values of u, 



6 ( Sl + E P( S2 + r 2 + r 3, S3 + r 3 ) J , 

s£ Z3„ = o \(r|u,k) J 

The sum over si is now trivial, so denoting t = (s 2 , S3), 
the expression simplifies to 

EE* E P( s 2 + r 2 + r 3 ,s 3 + r 3 ) 

teZ2 u=Q Y(r|u,k) 

The last step concerns the identification of the argument 
of $0 (the new measures for the effective system) . Using 
the translational invariant in t the above expression can 
always be written 

^^$ f Bl (k)p(t + e ; ) V (6) 

tez 2 "=o \/=o / 

where e = (0,0), ei = (1,0) and e 2 = (0,1), and the 
coefficients et/(k) are listed in Tabled (notice that $o(0) = 
0, so for some u and k there will be no contribution). 

It is now time to reinsert © back into Q, group terms, 
and write the final form for the functional, 



m\p]= E [*o(/>(s) + p(s + ei )) 



(7) 



+ * (p(s) + p(s + e 2 )) - 3$ (p(s)) 



The above expression can be more easily visualized by 
using the diagrammatic notation introduced in |17| . 



OTW = E [*o(c^) + $o(^)-3<I>o(< 

sez 2 



(8) 
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FIG. 2: Hard cubes with edge length a = 2 in a simple cubic 
lattice. Their centers of mass are constraint to lie on the plane 
si + S2 + S3 = 0. Filled circles mark the sites excluded by the 
cubes which lie on this plane. Notice that on the plane, cubes 
behave as the hard-hexagon lattice gas with nearest-neighbor 
exclusion. The diagram at the lower right corner helps to 
visualize the relative position of the cubes in the figure. 
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TABLE II: The same as Table Q for the linear combinations 
in 11 1 21 defining the weighted densities for the hard-hexagon 
lattice gas model. 



B. Triangular lattice 

As in the previous subsection, we will start with the 
FM functional for the parallel hard cubes with a = 2 in 
the simple cubic lattice. As the procedure we will follow 
is the same of the one used before, we will omit most 
details. Again, we will restrict the position of the centers 
of mass of the cubes to lie, in this case, in the plane 

V b = {( Sl ,s 2 ,s 3 ) eZ 3 : Sl + s 2 + s 3 = 0}. (9) 

A sketch of the effect of this confinement is shown in 
Fig. |21 It can be appreciated that the resulting effective 
lattice is a triangular lattice, and that the interaction 
becomes again nearest-neighbor exclusion. The effective 
system thus corresponds to the hard-hexagon lattice gas. 

The one-particle distribution function for this confined 
system can be written 



p(s%, s 2 ,s 3 ) = p(s 2 ,s 3 )5(si + s 2 + s 3 ), 



(10) 



where the arguments of the two-dimensional density pro- 
file correspond to the choice {ei = h 2 — bi, e 2 — h 3 — bi} 
for a vector basis. 

After eliminating the unnecessary degrees of freedom, 
the resulting contribution of (s) to the excess free- 



energy functional is 

E E $ ° ( E P(h + r 2 ,t2 + rs) J , (11) 

t6 Z 2 u=0 \(r|«,k) J 

using the same notation as in the previous case. Now 
u = r% + r 2 + r 3 . This can again be written 



E E*° (l>( k v( t + e ') 

tgz 2 «=0 



(12) 



J=0 



where the coefficients a;(k) appear in Table UTI (a third 
vector, ea = ea — ei, is involved this time). 

Inserting this in Q and regrouping terms, the excess 
free-energy functional for the hard-hexagon model turns 
out to be 



/^£M=E [$o(n (+) (s))+$o(n(->( 

sez 2 

-]T> (n«(s))+* (P(s)) 



(13) 



where we have defined the new set of weighted densities 
n (±) (s) = p(s) + p(s ± ei) + p(s ± e 2 ), 



nW(s) = p(s)+p(s + e i ), 1 = 1,2,3. 
In diagrammatic notation, 



(14) 



^] = Eh(A)+MV 

s6Z 2 

-$o(^)-*o(/) -*o(\) +*o(o) 



(15) 



C. Simple cubic lattice 

As in the two previous systems, we will write down a 
closed-form free-energy density functional for this model, 
but now we have to start from a four-dimensional system 
of a = 2 parallel hypercubes in a simple hypercubic lat- 
tice and constrain the centers of mass of the hypercubes 
to be in the three-dimensional hyperplane 

V sc = {(sx, s 2 , s 3 , Si ) G Z 4 : Sl + S2 +2s 3 +4s4 = 0}. (16) 

As this cannot be graphically sketched we should clarify 
why we choose this particular hyperplane. A hypercube 
with cr = 2 centered at the origin contains all lattice sites 
s with coordinates si = 0, ±1, I = 1, ... ,4. Now, the in- 
tersection of such a hypercube with "P sc can be split into 
two sets: one with S4 = and si + s 2 + 2S3 = 0, and the 
other one with s 4 = ±1 and si + s 2 + 2s 3 = =p4. These 
two sets define three parallel (three-dimensional) planes. 
The first one coincides with P S q, Eq. J3J| (hence the choice 
of the coefficients of these three coordinates in 'Psc) and 
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as Table Q] for the linear combina- 
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tions in 1231 defining the weighted densities for the nearest- 
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III: The same as Table for the linear combina- 
















tions in 11811 defining the 
neighbor exclusion lattice 


weighted densities for the nearest- 
gas in the sc lattice. 


where rS l \s) 


= p(b) + 


P(s + e z ), 


or diagrammatically, 



thus includes the five points (including the origin) of the 
two-dimensional square lattice of Sec. Ill Al The last two 
planes contain only one site each (symmetrically placed 
with respect to the origin). Excluding the origin, these 
sites complete the eight vertices of the octahedron defin- 
ing the 'shape' of the particle. 

The effective underlying lattice is expanded, e.g., by 
the vector basis of the hyperplane {ei = b 2 — bi,e2 = 
b 3 -bi-b 2 ,e 3 = b 4 -bi-b 2 -b 3 }, with {bi, b 2 , b 3 , b 4 } 
the canonical vector basis of Z 4 . Thus, the one-particle 
distribution function for this system can be written as 

p(s) = p(s 2 + s 3 + 2s 4 , s 3 + s 4 , 8i)S(si + s 2 + 2s 3 + 4s 4 ), 

(17) 

where, again, the three-dimensional density profile cor- 
responds to that of the effective system. 

Repeating the procedure described in the previous ex- 
amples we arrive at this expression for the contribution 
of n( k '(s) to the excess free-energy functional, 



££ $0 £ +M2 + r 3 + r 4 ,< 3 + r 4 ) 

tgz 3 "=0 \(r|u,k) 

where now u = r\ + r 2 + 2r 3 4- 4r 4 and v is a shorthand 
for r 2 + r 3 + 2r 4 . As in the previous cases, this becomes 



££$o I>(k)p(t + e«) , 



(18) 



tez 3 u=o \i=o 



with e = (0,0,0), ei = (1,0,0), e 2 = (0,1,0) and e 3 = 
(0, 0, 1), and the coefficients a;(k) given in Table ITTT1 

The final expression for the excess free-energy density 
functional of this model is 



rad = £ 



£<& (n«(s))-5$ (p(s)) 
.1=1 



, (19) 



mi\p} = £ [$ (^)+$o(c^) 

+0>o(c p )-5<i>o(o)]. (20) 



sSZ 3 



D. Face-centered cubic lattice 

The nearest-neighbor exclusion lattice gas in the fee 
lattice can be obtained from the same four-dimensional 
system we have used above, but now the centers of mass 
of the hypercubes are confined to the three-dimensional 
hyperplane 



{(si,s 2 ,s 3 ,s 4 ) G Z 4 : si+.s 2 +.s 3 + 2.s 4 = 0}. (21) 



A vector basis expanding the effective underlying three- 
dimensional fee lattice is {e 4 = b 2 — b 4 , e 2 = b 2 — b 3 , e 3 = 
b 4 — bi — b 3 }. When a a = 2 hypercube is placed 
at the origin (0,0,0,0) the excluded sites in V{ cc cor- 
respond to the set of twelve nearest neighbors of the fee 
lattice, whose coordinates in the chosen basis are {ei = 
(1,0,0), e 2 = (0,l,0),e 3 = (0,0,l),e 4 = (l,-l,0),e 5 = 
(1,0,— l),e6 = (0,1,-1)} and the opposite ones. Un- 
der this external potential the one-particle distribution 
function takes the form 

p(s) = p{s 2 + s 3 + s 4 , -s 3 - s 4 , s 4 )S{s 1 + s 2 +s 3 + 2s 4 ), 

where the density profile in the r.h.s. is expressed in the 
chosen basis. 

Proceeding as in the previous cases, the contribution 
of the weighted density n^(s) to the excess free-energy 
functional reads 

££ $0 f £ p( t i+ v M-r 3 -r 4 ,t 3 + r 4 ) J , 

tgZ 3 «=0 \(r|u,k) J 

where now u = r± + r 2 + r 3 + 2r 4 and v is a shorthand 
for r 2 + r 3 + r 4 . This is more conveniently rewritten 



££*" (£az(k)p(t + eOj, 

t e z 3 «=o \;=o / 



(23) 
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FIG. 3: Body-centered cubic lattice (black points) obtained 
by applying an infinite-strength external potential at the 
white points of a simple cubic lattice. Hard cubes with a — 2 
can only be placed at black points, thus the interaction po- 
tential becomes nearest-neighbor exclusion in the bcc lattice. 
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TABLE V: The same as Table[I]for the linear combinations in 
12811 defining the weighted densities for the nearest-neighbor 
exclusion lattice gas in the bcc lattice. 



sketched in Fig. |3| The effect of this external potential 
in the density profile of the original system amounts to 



with the coefficients a;(k) given in Table ITVl 

Gathering together the contributions of all weighted 
densities, as prescribed in , the excess free-energy den- 
sity functional of this model turns out to be 



E<*>o (n«(s))+5$ (p(s)) 



1=1 



(24) 



where 



n (±) (s) = p(s) + p(s±e 1 )+ p(s ± e 2 ) + p(s ± e 3 ), 
n {,) (s) = p(s) + p(s + ei), Z = l,...,6. 

(25) 

The diagrammatic notation is of little help in this case, 
so we omit it. 



E. Body-centered cubic lattice 

It is impossible to reduce a hypercubic four- 
dimensional lattice to a three-dimensional bcc lattice by 
projecting on a hyperplane. Therefore, the method used 
in the previous cases is not suitable for this one. But 
there exits a different procedure to obtain the functional 
for this system. Again, we will begin with a system of 
parallel hard cubes, but now of the same dimensional- 
ity of the target system. Also, the key ingredient of the 
method is to apply an infinite-strength external potential 
to the appropriate set of lattice sites. 

Let us consider a three-dimensional system of parallel 
hard cubes with a = 2 in a simple cubic lattice. If we 
now restrict the centers of mass of the cubes to the set 



M , s 2 + S 3 Si + S 3 Si + s 2 . , . , 
P(*) = P\ n ' o > n )Sc(s), (27) 



where #c(s) = 1 if s € C and otherwise, and the coor- 
dinates in the density profile of the r.h.s. are referred to 
the basis {ei = b 2 + b 3 h>i, e 2 = t>i + b 3 b 2 , e 3 = 
bi + b 2 b 3 }, {bi,b 2 ,b 3 } being the canonical vector 
basis in Z 3 . 

In order to calculate the contribution of each weighted 
density n( k )(s) we will take into account that 

E=E E - 

sGZ 3 u=0s+b„e£ 

where bo = (0, 0, 0). Then, after making use of the trans- 
lational invariance in s, we have for the contribution of 
n( k )(s), 



EE*o E P(s + r-b„)], 

sG£m=0 y(r|u,k) 

where condition (r\u, k) selects those r S B(k) such that 
r — h u G C. In a more convenient way, the previous 
expression can be written 



EE $ ° (E a <( k V( s + e <)V 

s6 Z 3 «=0 \l=0 I 



(28) 



where = ei+e 2 + e 3 , the coefficients a;(k) are given in 
TablelVl and all the spatial vectors appearing are referred 
to the vector basis of C. 

Taking into account the contribution of all the involved 
weighted densities, the total excess free-energy results 



£ = {(8i,8 2 , s 3 ) S Z 3 : all Si odd or all even}, (26) P^ c [p} = E 



E*o (n«(s))-7* (p(b)) 
.l—i 



(29) 



then C forms a bcc lattice and the exclusion of the cubes 
in the original lattice corresponds to nearest-neighbor ex- 
clusion in the effective lattice. This transformation is 



with n®(s) = p(s) + p(s + e/). Again, the diagrammatic 
notation is of no much help in this case either. 
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Lattice 


Order 


V 


Pp 


z 




square 


second 


1/2 


0.523 


1.687 


0.392 


triangular 


first 


0.684-0.754 


0.737 


7.70 




sc 


second 


1/3 


0.305 


0.763 


0.350 


fee 


first 


0.579-0.837 


0.461 


5.29 




bee 


second 


1/4 


0.216 


0.490 


0.305 



TABLE VI: Transitions of nearest-neighbor exclusion lattice 
predicted by FM functionals. 



III. THERMODYNAMICS 

In the previous section, we have obtained the density 
functional for all systems in closed form. From them it 
is possible to derive all the equilibrium properties of the 
system in the presence of an arbitrary external poten- 
tial. In this section, we will restrict ourselves to the bulk 
phase diagram. The description of bulk behavior that 
lattice FM functionals provide turns out to be equiva- 
lent to that obtained from other classical theories. Thus 
in this particular application it provides nothing new, 
and we include it both for completeness and to show the 
remarkable fact that FM theory subsumes many other 
theories in a unified framework. Anyhow, we want to 
stress that what FM theory provides are functionals, and 
functionals perform best when applied to inhomogeneous 
problems. 

The phase transitions obtained for each model are col- 
lected in Table EH 



A. Square lattice 

As the interaction is only between nearest neighbors, 
we will distinguish two sublattices in such a way that the 
nearest neighbors of every site of one sublattice belongs 
to the other sublattice. Let the density of each sublattice 
be p\ and p 2 - Given a total density p it holds 2p = p\+p2- 
When we consider such a density profile in J7J the free- 
energy density becomes 

$ sq = $id + 2$ (?7) ~ |[$o(pi) + $o(r) ~ Pi)], (30) 

where < rj = 2p < 1 is the packing fraction, = 
(1/2) J2 . pi (In pi — 1) is the ideal contribution and we 
have written p 2 = rj — p\ . 

The equilibrium state of the system at constant pack- 
ing fraction r\ is given by the global minimum of l|3UI) with 
respect to p\. This is a solution to the Euler-Lagrange 
equation 



Pi(l-Pi) 3 = {rj - Pl ){l - r, + Pl f 



(31) 



For low values of r\ this state corresponds to the uni- 
form or disordered one, with p\ = pi = P- For higher 




FIG. 4: Phase diagram (pressure vs. packing fraction) of the 
hard-square lattice gas in the lattice FM approximation (solid 
line). The system undergoes a second-order phase transition 
from a disordered state to an ordered one at rf^ = 1/2. Dot- 
ted line: results of the 'ring' approximation obtained by Bur- 
ley m. 



values, one expects the system to undergo a phase tran- 
sition to an ordered phase, where one of the sublat- 
tices is preferentially occupied. This situation is ex- 
actly the one described by the solutions of Eq. I|31[) . For 
< rj < r/l q = 1/2 the minimum is given by p\ = p. 
Equation of state and fugacity of the disordered phase 
can be easily obtained from i|3U[l . 



hi 



(1-W2) 3 



(i-v) 2 ' 

For 1/2 < rj < 1 , we have 



'fluid 



?7(l-r7/2) 3 
2(1 ~r,Y 



eq x 

P "=2 



V + (2 - V) 



277 —1 
3-2?/ 



(32) 



(33) 



This gives the ordered state. Equation of state and fu- 
gacity in this phase are given by 



solid 



3- 2rj 
4(1 - t?)' 



Pi q (l-Pi q ) 3 

(i - v y 



sq 

2 solid — 



(34) 

The equation of state of both phases is plotted in 
Fig. Also, the occupancy of each sublattice is plot- 
ted in Fig. The transition is second order, and the 
critical point values are listed in Tabic IVTl These are the 
same as those obtained by Burley [23j and Temperley 
|24| with a Bethe approximation. The next correction 
to Bethe approximation (called 'ring' approximation ) in 
the scheme proposed by Rushbrooke and Scoins |25ll2l| 
has also been considered by Burley j2^|. The results are 
plotted in Fig. 0] for comparison. The qualitative phase 
behavior is the same, but the critical values improve. A 
further correction ('necklace' approximation) within the 
same scheme has been implemented by Temperley |27j . 
who claims that the transition becomes first order. 
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aT 0.5- 





FIG. 5: Densities of the two sublattices, pi and p2, for the 
hard-square lattice gas, as a function of the packing fraction 



B. Triangular lattice 

As a triangular lattice has three sublattices such that 
a nearest neighbor of a site in one sublattice belongs to 
the another sublattice, we will characterize the density 
profile with the densities of each sublattice: p\, p2 and 
P3, the total density of the system being related to them 
through 3p = pi + p 2 + P3- 

The free-energy density for a three sublattice configu- 
ration can be obtained from 

$ tr = *id + 2$ (??) - *o(pi + pi) - $o(pi + Pi) 
- $o(p 2 + Pa) + g NOi) + *o(P2) + $oGo 3 )] , (35) 

where < r\ = 3p < 1 and $id = (1/3) J2i Pi( m Pi — !)■ 
As in the previous case, the phase diagram can be ob- 
tained by minimizing this free-energy density at constant 
packing fraction. For low values of the density the stable 
phase is uniform, so pi = P2 = P3 = P- In contrast, for 
high values of the density one sublattice is preferentially 
occupied. In this ordered phase the configuration that 
results from i|35|) is p\ > p 2 — ps- Thus, in what follows 
we will consider the latter two sublattices as equivalent. 
The Euler-Lagrange equation obtained from (|35|) is 



p 2 (l- 2p 2 ) 3 (l ~y + 2p 2 ) 

(l-2p 2 )(T ? -2p 2 )(l-7 7 + p 2 ) S 



(36) 



The solutions to this equation as a function of r\ indicate a 
first-order transition. The coexisting packing fractions as 
well as the pressure and fugacity at coexistence appear in 
Table IVTl The equation of state for the disordered phase 



fluid 



hi 



(1 - 2,y/3) s 



(l-r / /3)(l- ?7 ) 2 ' 



(37) 



FIG. 6: The same as Fig. 2] f° r the hard-hexagon model 
(nearest-neighbor exclusion in a triangular lattice). Solid line 
represents the equation of state from the FM approximation. 
Dashed line represents the exact solution by Baxter (2S)l . In 
our approach the system undergoes a first-order phase tran- 
sition from a disordered phase to an ordered one. Coexisting 
packing fractions are ?7fl u id = 0.684 and r; so iid = 0.754. It is 
worth mentioning the good agreement of FM approximation 
with the exact solution at low and high densities, despite its 
failure in the critical region. 



and the fugacity 



7/(1 - 2?7/3) e 



" fluid 3(1-t?/3)(1-t7) 6 " 
For the ordered phase 



and 



"solid 



(l-?? + p 2 ) 6 (l-2p 2 ) 3 

(l- ?? )6(l_^ + 2p 2 )(l-p 2 )2' 



P2{i-v + P2?{i-2p 2 y 

(1-^) 6 (1-P2) 



(38) 



(39) 



(40) 



where p 2 is given by l|3ofl . The phase diagram is plotted 
in Fig. |SJ Sublattice densities are shown in Fig. [7| 

This result coincides with that obtained by Burley [2j| 
by using the 'ring' approximation (see Sec. IIII A|l . and 
was also obtained by Tcmperlcy through a combinatorial 
finite matrix method [24| . 

This system has been exactly solved by Baxter p^.l3fj| 
(the exact solution is shown in Figs. El and [7J. The 
system exhibits a continuous phase transition (the crit- 
ical point in the pressure vs. density diagram is a hor- 
izontal inflection point) at packing fraction 7y° xact = 
(3/10)(5 - V5) « 0.829 and critical pressure f3p c c xact = 
(l/2)ln[(27/250)(25 + ll%/5)] w 0.839. This is the first 
time that FM theory can be compared with an exact re- 
sult in dimension higher than one. Although the theory 
fails in predicting correctly the order of the transition, 
the agreement in the whole range of density except the 
critical region is rather accurate. 
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FIG. 7: The same as Fig.|5]for the triangular lattice. Full cir- 
cles represent the configuration of coexisting phases. Dashed 
lines are sublattice densities in the exact solution. Empty 
circle is the exact critical point. 



FIG. 8: The same as Fig. [I] for the nearest-neighbor exclu- 
sion lattice gas in the simple cubic lattice. A second-order 
transition is predicted at r)* c = 1/3. 



C. Simple cubic lattice 

In the close-packed state the particles of this system 
occupy one sublattice with fee symmetry, while the sub- 
lattice formed by nearest neighbors of that one remains 
empty. Thus, we will only consider these two sublattices, 
with densities p\ and p2, respectively. Again, the total 
density p is related to these numbers by 2p — p\ + p2 ■ 

When we insert such a density profile in the functional 
(|19fl we get for the free-energy density 



$ sc = $ id + $o(2p) 



2 I $ o(Pi) 



$o(2p- Pl )], (41) 



where we have eliminated the dependency on p^. For a 
fixed value of the density, the global minimum of i|41|) is 
a solution to the Euler-Lagrange equation 



iv - pi)(i - v + pi)' 



= i, 



(42) 



< r\ = 2p < 1 being the packing fraction. When < 
V < = 1/3 j the minimum is given by the uniform 
phase = r\. For the equation of state and fugacity we 
have 



/Muid = In 



(1 - rj) 3 



_ 77(1 - 77/2)5 
Zfluid " 2(1-77) ■ { ' 



Upon increasing the density the system experiences a 
second-order disorder-order transition at r/j? c = 1/3. For 
1/3 < 77 < 1, the density of the preferred sublattice can 
be calculated from l|42|l to be 



eq 

Pi = 



77 +(2 -77)1 



'2[5 - 477(2 - 77)] 1 / 2 - 5(1 - 77) 



377 



(44) 



or 0.5 




FIG. 9: The same as Fig.|K|for the simple cubic lattice. 



In this phase, the equation of state is given by 



solid 



(3d 



fluid 



hi 



2(5 - 477) - y/5 - 477(2 - 77) 



5-377 



and the fugacity by 



Pi q (l-Pi q ) 5 
(1 - 77)6 



v solid 



(45) 



(46) 



The equation of state is plotted in Fig. |SJ Sublattice 
densities are shown in Fig. The critical point values 
are listed in Table IVll As for the square lattice, for this 
model the result from FM theory coincide with a Bethe 
approximation 23] . 
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FIG. 10: The same as Fig. 0] for the nearest-neighbor exclu- FIG. 11: The same as Fig. |S]fbr the fee lattice. Solid circles 

sion in the face-centered cubic lattice. A first-order phase represent configurations of coexisting phases. The inset shows 

transition takes place at pressure f3p = 0.461 with coexisting a detail of sublattice density p2. 
states at rj^ id = 0.579 and rf™ m = 0.837. 



and 



D. Face-centered cubic lattice 

An fee lattice contains four equivalent sublattices each 
with a sc symmetry. As in the previous cases sublat- 
tice densities are denoted, pi,p2,P3,Pi, and rj = Ap = 
Pi + P2 + P3 + Pi holds for the total density p or packing 
fraction 77. 

The free-energy density for such a density profile in the 
FM approximation is obtained from 1241) as 

e 4 

i<j »=1 

where <I>id = (1/4) Y^i PiQ- n Pi — !)■ The global minimum 
of $f cc at fixed density with respect to the sublattice 
densities yields the thermodynamically stable phase. But 
now there is a subtle point we have to take into account: 
the close-packing is degenerated. If we consider the fee 
lattice as a stacking of square lattices each having its 
sites on the centers of the squares of the previous one, 
then we can fill alternative layers independently. Thus, 
we will consider different densities for the four sublattices 
and the minimum of the free energy will determine the 
structure of the most stable phase. (Previous studies 
|3lj impose the equivalence of three sublattices from the 
beginning, i.e. pi ^ p 2 = pa = Pi-) 

Solving numerically the Euler-Lagrange equations ob- 
tained from (|47|) at fixed packing fraction, we obtain the 
equation of state plotted in Fig. ^3 Sublattice densities 
at equilibrium phases are shown in Fig. 1111 The system 
undergoes a first-order phase transition to an ordered 
phase with p\ ^ p 2 = pz = p±. (See Table fVTI for the 
critical point values.) For the uniform phase 



In 



(1_„/4)5(1_„)2 



(48) 



_fcc 



fluid 4(l-77/4)5(l-7?)8' 
For the ordered phase 

{&» -V (l-^ + 2p 2 ) 12 (l-2p 2 ) 12 

PPS0M 4 (l-7 ? + 3^)5(1-^)15(1-77)8 

and 



fee _ p 2 (l-r/ + 2p 2 ) 4 (l-2p 2 ) t 

z solid — 



(l- P2 )5(l-77)8 



(49) 



(50) 



(51) 



with p2 the minimum of the free energy 147J1 . 

It is worth mentioning that a metastable phase transi- 
tion from the disordered state to a smectic ordered state 
(alternative square-lattice layers are uniformly occupied, 
with densities p\ and p 2 ) is also obtained from the func- 
tional P7|) . 



E. Body-centered cubic 

At close packing state this system has one sublattice 
completely filled while the nearest-neighbor sublattice is 
empty. Therefore, we will consider the density profile as 
in the previous similar cases (square or sc lattices). 

The free-energy density [Eq. I|29|)] becomes 

$bcc = $. d + 4$o(?7) _ 7 ^ o{pi) + <Pa{v _ pi)] . (52) 



The Euler-Lagrange equation is 
Pi(l-Pi) 7 



iv - Pi)(l - V + Pi)"' 



(53) 
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FIG. 12: The same as Fig.[I]for the nearest-neighbor exclusion 
lattice gas in the bcc lattice. The system undergoes a second- 



order phase ordering transition at 77° 



1/4. 



The phase behavior obtained from this equation is sim- 
ilar to that of the square or the sc lattices: The system 
undergoes a second-order phase transition at ?yj? cc = 1/4; 
below this density the stable phase is a fluid, while above 
the particles occupy preferentially one sublattice. The or- 
dered solution of ijSlfl) can be obtained analytically, but 
the expression is rather cumbersome, so we just plot it in 
FigED For the fluid phase, equation of state and fugacity 
are given by 



Muid = In 



(l-^/2) 7 



(1 - r,)4 

For the ordered phase, we have 



bcc 
z fluid 



77(1 - 77/2)' 

2(1-77)8 



ftPsolid 



— m ; r-3 . 

2 (l-^) 4 



and 



bcc 
"solid 



Pl(l-Pl) 7 
(1-77)8 > 



(54) 



(55) 



(56) 



with pi the ordered solution of JHSJ). 

These results are again equivalent to those obtained 
with a Bethe approximation [23l | . The equation of state 
is plotted in Fig. ^3 (See Table IVT1 for the critical point 
values) 



IV. DISCUSSIONS AND CONCLUSIONS 

Dimensional crossover is a property that connects dif- 
ferent models. A system of parallel hard cubes (both on 
and off a lattice) constrained to lie on a plane parallel 
to their sides becomes a system of parallel hard squares. 
A system of parallel hard squares constrained to lie on 
a straight line parallel to their edges becomes a system 



FIG. 13: The same as Fig.|S]for the bcc lattice. 



of hard rods. These are obvious dimensional crossovers, 
but in this paper we have introduced a few more: paral- 
lel hard cubes constrained to lie on more general planes 
produce nearest-neighbor exclusion lattice gases in the 
square and triangular lattices; four-dimensional paral- 
lel hard hypercubes constrained to lie on certain hyper- 
planes become nearest-neighbor exclusion lattice gases 
in the sc or fee lattices. If this dimensional constraint is 
regarded as the application of an infinite-strength exter- 
nal field, then more general patterns of such an external 
field produce new models out of the parallel hard cube 
one (the nearest-neighbor exclusion lattice gas in the bcc 
lattice, for instance). 

The message we want to transmit with this work is 
that FM functionals for all these models have the remark- 
able property of being connected with each other through 
the same transformations. More than that: further di- 
mensional crossovers of these models are also consistently 
captured by FM functionals. For instance, all functionals 
obtained in this paper produce the exact functional for 
hard rods when reduced to one dimension; also, the func- 
tional for the nearest-neighbor lattice gas in the square 
lattice, Eq. J7J, can be obtained from the one in either 
the sc lattice, Eq. (|T§j) . the fee lattice, Eq. or the 

bcc lattice, Eq. (|29(l . through the application of an ap- 
propriate external potential. Similarly, the functional for 
the nearest-neighbor lattice gas in the triangular lattice, 
Eq. (|13fl . can be recovered from that in the fee lattice, 
Eq. I|24|) . Such a degree of internal consistency is not 
shared by any other known density functional theory, and 
it puts FM theory at a different level. As a byproduct, 
as pointed out in the Introduction, it warrants a good 
behavior of FM functionals when dealing with highly in- 
homogeneous situations. 

When applied to study the bulk phase behavior of 
these nearest-neighbor lattice gases, FM functionals pro- 
duce reasonable results, sometimes with important dis- 
crepancies in the critical region, sometimes more accu- 
rate, and always very accurate at low and high densities. 
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There are, no doubt, better methods to fit the equation of 
state and obtain a better description of this phase behav- 
ior, like, for instance, finite-size analysis or 
series expansions |3l|, yj, yg, yTj • They are particularly 
useful in predicting the critical behavior. FM approach 
has an important advantage on these method, whatever 
the sacrifice in accuracy: it leads to simple, closed-form 
functionals. Thus it permits to study inhomogeneous 
problems, which are absolutely out of the scope of these 
other more accurate methods. A particularly interesting 
example of such inhomogeneous problems can be found 
in recent studies of fluids in porous media 0, l^jj, 0] , 
for which lattice models seem to capture enough physical 
information to describe many interesting phenomena not 
yet described by other models. 

An interesting observation to make from the results of 
the application of FM theory to describe bulk phase be- 
havior is the close connection it has with other classical 
approaches. For those lattices with only two sublattices 
(loose-packed lattices) FM theory reduces to a Bethe ap- 
proximation (square, sc and bcc lattices); for lattices with 
more than one sublattice (close-packed lattices) it be- 
comes equivalent to another cluster-like approximation. 
In fact, when closely looked at, all these approximation 
are but particular cases of Kikuchi's cluster variation the- 
ory |S8j . This theory proposes a hierarchical scheme of 
approximations under the basis of describing in an exact 



manner clusters of increasing size. Roughly speaking, 
the larger the cluster the more accurate the results — and 
the more involved the theory. The connection between 
FM theory and the cluster variation method is definitely 
worth exploring, because one drawback of the latter is 
that there is not an a priori criterion to choose the clus- 
ters, other than "the larger the better"; as a matter of 
fact, some clusters (depending on the model) produce an 
optimal result and some others (even larger ones) spoil 
the accuracy, and the reason for that is unknown. FM 
theory, on the contrary, leaves no freedom to choose the 
clusters, but those it prescribes seem to be optimal in the 
above sense. It is very illustrative Burley's treatment of 
the nearest-neighbor exclusion lattice gas in the fee lat- 
tice |23j : the inappropriate election of the clusters [only 
part of those prescribed by the FM functional J23J] leads 
to a spurious divergence of the free-energy at a density 
lower than the close-packing. This is a line of investiga- 
tion we are currently following. 
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